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Abstract 

This paper gives suggested values for some of the input parameters used in the CAP-TSD 
computer code. These parameters include those associated with the mesh design and time step. 
The guidelines are based principally on experience with a one-dimensional model problem used 
to study wave propagation in the vertical direction. 

Introduction 

Wave propagation through the finite-difference meshes employed in CFD may be simulated 
most simply by one-dimensional model problems. For TSD theory, an appropriate model is the 
second order wave equation, which describes the vertical waves in the CAP-TSD code (ref. 1). 
The use of this equation provides significant advantages: the exact solution of the partial dif- 
ferential equation is well known for any initial/boundary conditions; the numerical solution of the 
finite-difference version of the continuous equation is very inexpensive to compute; analysis of the 
difference equations may be used to aid in interpreting the numerical results; and computer graph- 
ics animation of the wave motion is easy to produce and clearly demonstrates wave propagation, 
reflection, and distortion effects. 

In the transonic small disturbance code CAP-TSD independent one-dimensional meshes are 
used in each of the computational coordinates £, rj, and (. Typically, the streamwise £ and 
vertical ( meshes extend 10 to 40 chord lengths from the configuration and the spanwise rj 
mesh extends 2 to 4 span lengths. The choice of number, extent, and spacing of mesh points is a 
difficult problem. Once these choices have been made, the interactive computer code described 
in reference 2 may be used to aid the mesh designer. 

Results obtained with the simple model problem were reported in reference 3. The algorithm 
employed was that for the vertical sweep in the CAP-TSD code. Wave reflection and time step 
effects similar to those experienced with two-dimensional TSD codes were observed. Analysis 
of the finite-difference equations supported the numerical results. Additional calculations for this 
model problem are reported herein. 


The guidelines presented in this paper are based on three levels of experience: extensive 
calculation and analysis for the simple model problem; use of the CAP-TSD code, especially for 
two-dimensional flow; and physically plausible requirements for modelling the transonic flow All 
calculations and figures were produced using a SUN workstation. 


Symbols 


h 

^far— field 

^m«* 

^surface 

k 

U 

'''max 

M 

N cy 

N p C 
Np W 

N z 

" max 

P 

t 

tc 

T 


z 

^mu 

P 

At 

A 


local mesh spacing 

mesh spacing at the far-field boundary 
maximum mesh spacing 
mesh spacing near the wing surface 
reduced frequency 

maximum reduced frequency for the mesh 
Mach number 

number of cycles of harmonic motion 
number of time steps per cycle 
number of mesh points per wave length 
total number of mesh points 
velocity potential 
time 

time at pulse center 
temporal period 

total time for converged calculation 
streamwise coordinate 
vertical coordinate 

maximum coordinate (far-field boundary location) 

VI - M 2 

time step size 

spatial wave length 

computational coordinates 


All quantities are normalized — lengths by chord, time by chord over stream speed, and frequency 
by stream speed over semichord. 


Wave Propagation 

In two-dimensional linear theory the wave fronts are circles which propagate outward with 
sonic speed as their centers are convected downstream at the stream speed. These wave patterns 
for four subsonic Mach numbers are shown in figure 1 . In each case the mesh extends 20 chord 
lengths from the origin both vertically and streamwise. The circles represent the wave crests 
created by a harmonically oscillating disturbance located at the origin. Several comments on 
these results will be helpful in understanding the mesh and time step sizes needed to describe 
wave propagation using a TSD code. 
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Since time is normalized by stream speed, in one cycle the centers of the circles are convected 
downstream a distance that is independent of the Mach number. Attention is drawn to the smallest 
circle in each figure to illustrate this conclusion. The center of this circle lies at 

X —T = Tr/k = 2 7T 


in which the first equality results from the normalization of t by stream speed and chord, the 
second from normalization of k by semichord, and the third from the choice of k = 0.5 for this 
case. 



Figure 1- Wave patterns for two-dimensional linear theory. 
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The radii of the circles depend on the Mach number. The wave length as measured from the 
center of the circle is 

A = T/M — n/Mk 

which decreases with increasing Mach number. With time normalized by stream speed, waves 
propagate more quickly through the mesh as the Mach number is decreased. Therefore, for a 
fixed total calculation time T max , disturbances travel a greater distance, which leads to a greater 
opportunity for reflections from the boundary to affect the solution. 

On the other hand, as the Mach number is increased the wave pattern as viewed from the 
origin becomes more and more distorted. Waves propagate slowly upstream and the effective 
wave lengths that must be resolved on the mesh become shorter, implying a need for very fine 
mesh spacing to track these disturbances accurately, or conversely, leading to increasingly large 

errors on a fixed mesh due to inadequate resolution. Specifically, wave lengths for disturbances 
propagating along the axes are 


1 -Mir 

nr i upstream 

1 + M 7C 

— — — — downstream 


X= Mk VertiCal,y 

Since propagation speed is directly proportional to wave length, the time required for a disturbance 
to reach a far-field boundary will be inversely proportional to the wave length. 

The remainder of the paper contains three main sections. The first section presents results 
or a pulse disturbance which illustrate the effects of time step and mesh spacing on frequency 
response. The next section contains the recommendations for time step and mesh properties. The 
final section shows results for harmonic disturbances which illustrate the efficacy of the guidelines. 

Pulse Results 

All pulse calculations were made on the same mesh. This mesh has a spacing which is 
fine near the origin, stretches smoothly to a maximum, and then shrinks somewhat at the far-field 
boundary. The mesh is defined by z as a fifth degree polynomial function of the mesh index and 
was designed using the program described in reference 2. The mesh properties are 


^"max — 20 

= 32 

^•urface — 0.01 

^far-field = 0.5 
flmaj = 1.0123 

where the mesh spacing h would equal A z for the vertical mesh in TSD theory. 

A Gaussian pulse in angle of attack is routinely used to produce the frequency response 
function in applications of the TSD code. Such a pulse contains energy at all frequencies; however, 
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most of the excitation occurs at low frequency and the effective range depends on the pulse width. 
The pulse responses for four time step sizes are presented in this section. The calculations give 
the solution for the wave equation with unit speed 


dt 2 

for an input in downwash centered at t e given by 



applied at z = 0. The function P represents the potential in TSD theory. The exact solution is 

It may be remarked that for a vertical wave the unit propagation speed assumed here occurs for 

M = p = Vol ~ 0.7 


Figure 2 gives the results for a large time step size. The lines with symbols are the numerical 
results and the plain lines are the exact solution. The discussion of the results for this figure will 
be given after some general comments on its organization. 

The lower part of the figure shows the time history of the solution at z = 0. In the TSD 
code this solution would correspond to the pulse in potential (or pressure, or lift) which results 
from the input pulse in angle of attack. The heavy tic-mark near t = 42 marks the time for the 
initial disturbance, centered at t c = 2.36 (the third time step for this case), to travel to the far-field 
boundary and return. 

The figure in the upper right gives the shape of the wave in space at a time shortly before 
the propagating pulse reaches the far-field boundary located at Z mftX = 20. The tic-marks show 
the mesh points, with fine spacing at z = 0, moderate spacing at z = 20, and coarse spacing 
near z = 10. 

The figure In the upper left gives the Fourier transform of the time history shown at the bottom. 
The real and imaginary parts of the transform are plotted against reduced frequency for k up to 
2.0. Most of the energy for a pulse of the width used is contained in this frequency range. 

Details of the parameters for this case are printed in the upper right of the figure. In this list 
dt is used for A t and CN is the ratio of time step size to mesh spacing. 
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Figure 2 gives the pulse response for A t = zr /4. This time step is too large, with fewer than 
four steps contained in the pulse, as may be seen in the time history at the bottom. The resulting 
severe broadening and loss in amplitude of the propagating pulse are seen at the upper right. 
This dissipation error is a consequence of the large time step. The frequency response function 
at the upper left is very smooth, but is accurate only at the lowest frequencies. Without the exact 
solution for reference the errors at higher frequency would not be detected. 


Implicit, 2nd order dt 
dt = 0.78540 tc = 2.36 
Newton = 1 

0.0100 < dz < 1.0123 

0.7758 < CN < 78.5398 

sig = 2.0 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 80 



MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, dt = 0.785 


Figure 3 gives the pulse response for Af = ir/8, half of the value used in the preceding 
figure. The frequency response results are fairly good up to about k = 0.4, but the propagating 
wave still suffers from dissipation errors. Dispersion errors, the spreading out of the frequency 
components as the wave propagates through the mesh, are also present. As will be apparent 
shortly, the latter errors become more severe as the time step is reduced, that is, as the spatially 
discretized equations are solved more accurately. 


Implicit. 2nd order dt 
dt = 0.39270 tc = 2.75 
Newton = 1 

0.0100 < dz < 1.0123 
0.3879 < CN < 39.2699 
sig = 2.0 

Zmox = 20. Nz = 32 

Tmax = 63. Nmax = 160 



MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, dt = 0.393 


Figure 3- Pulse response for At = 0.393. 
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Results for a further reduction in time step to At = tt/ 16 are given in figure 4. The dissipation 
errors in the upper right part of the figure are smaller but the dispersion errors are larger as 
compared with those in figure 3. These errors result from too few mesh points to represent 
adequately all of the wave lengths present in the solution and lead to shorter waves (high frequency 
components) propagating more slowly than longer waves (low frequencies). This spreading out 
of frequencies as the wave propagates is a consequence of the spatial discretization inherent in 
the finite-difference method and limits the frequency resolution attainable. The shortest waves 
that may be described on the mesh have two mesh points per period for which A = 2/Jfc. These 
dispersion errors also affect both the time history and the frequency response. There are only 
about three points contained within the width of the pulse in the coarsest part of the mesh. 

Oscillations in the time history occur as the wave propagates into the region of increased 
mesh spacing and consequently a decreased number of points per wave length. These errors 
reflect back to the left and contaminate the solution at the origin. 


MODEL EQUATION FOR CAP-TSD Implicit. 2nd order dt 

dt = 0.19635 tc = 2.95 

Quintic 32 pts, dt = 0.196 Newton = 1 

0.0100 < dz < 1.0123 
0.1940 < CN < 19.6350 
sig = 2.0 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 320 



The frequency response is reasonably accurate up to about Jfc = 1.0, although some small 
errors may be detected earlier. The catastropic breakdown in the solution above k = 1.0 arises 
because this mesh has a maximum spacing = 1.0123 for which two points per wave length 
implies A = 2.0246. Consequently, the highest frequency which may be represented is 

^max = 2/A ~ 1.0 

This value fc max is indicated with a heavy tic-mark on the figure. 

Figure 5 gives the pulse response for At = tt/ 32, the value recommended in the guidelines 
below. The time history of the response at the origin (airfoil) shows the severe oscillations that 
result from trying to propagate frequencies that are too high for the maximum mesh spacing used. 
However, these large errors do not affect adversely the low frequency response. Although not 
shown here, further reduction in At leads to little change in the quality of the solution. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, dt = 0.098 



Implicit, 2nd order dt 

dt = 0.09817 tc = 2.95 

Newton = 1 

0.0100 < dz < 1.0123 

0.0970 < CN < 9.8175 

sig - 2.0 

Zmax = 20. Nz = 32 

Tmax - 63. Nmax = 640 



0 z 20 




Figure 5- Pulse response for At = 0.098. 
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Guidelines 


Guidelines for mesh and time step parameters, with an indication of the reasons for their 
selection, are presented in this section. These choices are motivated by the results of the previous 
sections. Although the study leading to these recommendations strictly applies only to the vertical 
mesh, the guidelines should be reasonable for the streamwise mesh as well. As the transonic Mach 
number approaches one, the shrinking wave lengths illustrated in figure 1 for both the upstream 
and vertical waves would suggest that no mesh could be fine enough to properly treat any but the 
lowest frequencies. In the guidelines which follow these features are not addressed, and therefore 
the recommendations strictly apply only to the vertical mesh at M « 0.7. 

Mesh parameters 

Choice of the extent and spacing of the vertical £-mesh is based on the following consider- 
ations: 

(a) the far-fiekf boundary should lie at least 20 chord lengths from the wing; 

(b) to provide accuracy in both the near-field solution and the downwash boundary condition, the 
first mesh point should lie at about 0.00001 and the mesh spacing at the wing surface should 
be comparable to that used chordwise of 0.002 to .02 (for unit chord); 

(c) to avoid boundary reflections, the spacing at the far-field boundary should allow at least 8 
mesh points per wave length at the maximum frequency of interest; and 

(d) to avoid internal mesh reflections (i.e., to reduce dispersion errors), the maximum spacing 
should provide at least 4 points per wave length at the maximum frequency. These consid- 
erations lead to the following constraints on the mesh spacing h : 

^surface ~ 0.01 
^■far— field = 7r/8 k 

^mai = 7r/4A; 

in which A iurface is the spacing between the first two mesh points at the surface. The last 
two requirements result from the approximate relation (for a wave of unit speed) between the 
spatial wave length A and the number of mesh points per wave length N pw 

X = hN pv/ ss n/k 


from which 


h « 7r / 7V pw k 


Time parameters 

Selection of the time step size A t (DT in CAP-TSD) and number of time steps N mux 
(NSTEP) is based on the following considerations. The time step must be small enough to resolve 
the frequencies of interest and to capture significant transients. The total time T max must be great 
enough to allow for transients to be convected away from the configuration and for the steady-state 
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solution to develop. For an oscillatory case, at least two cycles of motion N cy should be used 
with at least 64 time steps per cycle N pc . For a reduced frequency k the period is 


T = iv/k 


The time parameters are selected such that 


At < 0.1 
r max > 20 tt 

N cy >2 

N pc > 64 

For a given frequency k the following algorithm may be used to meet these criteria: 


1. 

N p C = 

32/Jfc 



2. 

if N pc 

< 64 

then 

N pc = 64 

3. 

U 

a? 

'tT 
1 ll 

< 

4. 

N cy = 

20k 



5. 

if N cy 

< 2 then 

N cy =2 

6. 

N 

1 f max 

= N cy 

N pc 


7. 

T 

J max 

= N ma , x At 



This algorithm produces the following input values for different frequency ranges 

For 0.0 < k < 0.1 At = tt/ 32 and AT max = 64/A: 

For 0.1 < Jfc < 0.5 At = tt/32 and JV max = 640 

For 0.5 < k At = 7r/64 k and JV max = 1280A: 

It is noteworthy that JV max attains its minimum of 640 in the frequency range of greatest interest. 


Harmonic Oscillation Results 


The results presented here for the model problem with harmonic input illustrate the effective- 
ness of the guidelines given above. Calculations are first presented for a range of frequencies 
using the mesh described previously. According to the guidelines, this mesh should be adequate 
for frequencies up to about k — 0.75. Finally, a new mesh is designed for use at the maximum 
frequency of Jfc = 1.5 to verify the adequacy of the recommendations. 

A sequence of eight values of frequency is presented for k ranging from 0.025 to 1.5. The 
calculations give the solution for the wave equation with unit speed (M = /? « 0.7) 



d 2 

dz 2 


P 


for an input in downwash given by 


— P(0,f) = 2ksm2kt 
(J z 
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applied at z = 0. The exact solution is 


P(z,i) = cos2Jfe(f — z) 

The layout of the remaining figures is very similar to that used earlier, with two exceptions 
(refer to figure 6). Now, the Fourier components are computed from the last cycle of the time 
history shown in the bottom figure. The mean and first four harmonics are plotted with symbols 
and are listed in the table. For an exact result, all entries would be 0.0000 except for a 1.0000 in 
the real part at the input frequency (k = 0.025 for this case). The exact transform shown is the 
continuous result obtained by integrating the cosine response over one cycle. Also, the picture of 
the wave as it propagates across the mesh at the upper right is now shown for t = T max . 

We now turn to a discussion of the results in the order of increasing frequency. In each case 
the time parameters were chosen using the guidelines given above. The mesh is held fixed until 
the final figure. 
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The results for k = 0.025 are shown in figure 6. As is typical for low frequencies, the time 
step is at its upper limit, which requires a large number of steps per cycle and, consequently, 
a large number of time steps to obtain two cycles of oscillation. These large numbers are not 
required for accuracy in the model problem, but are appropriate for a TSD calculation. Notice 
that only a small portion of one wave is present across the spatial mesh, and the mesh is much 
finer than is required to represent the wave accurately. This mesh fineness leads to very small 
truncation errors and the Fourier components are essentially exact. The small differences shown 
may be round-off errors due to the 32-bit word length of the SUN workstation used. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, k = 0.025 

k ReF ImF 

0.000 -0.0006 0.0000 
0.025 1.0003 0.0001 



Implicit, 2nd order dt 

dt = 0.09817 k = 0.025 

Newton = 1 

0.0100 < dz < 1.0123 

0.0970 < CN < 9.8175 

Npc = 1280 

Zmox = 20. Nz = 32 

Tmax = 251. Nmax = 2560 




Lilli i i i i i i i i i i i i i i i i i 1 1 1 nl 
0 z 20 



Lilli mu mi Jiniiiiiiii iiiiiiiiiiin Ill mil ninnnJii n mil ulllU 11111 iiiiiiii III! Illl ill lllll llll III! Mill Illlllllllll lllll ill llllllllli 

0 t 252 


Figure 6- Harmonic response for k = 0.025. 
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The results for k = 0.05 are shown in figure 7. The quality of the calculation is again quite 
high. For this case the number of steps per cycle is N pc = 640. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, k = 0.050 


Implicit, 2nd order dt 

dt = 0.09817 k = 0.050 

Newton = 1 

0.0100 < dz < 1.0123 


k ReF ImF 

0.000 0.0003 0.0000 

0.050 1.0015 - 0.0004 

0.100 0.0000 0.0000 



0.0970 < CN < 9.8175 
Npc = 640 

Zmax = 20. Nz = 32 

Tmax = 126. Nmax = 1280 




Lmintiimiiimi mil 

0 


I II 1 1 II 11 1 ,1111 1 1 I I 1 1 1 HII II 1 1 II n n n i ill i n n ii ii ii mi n n i i ii H 1 1 in || || nil 1 1 1 m 1 1 


t 


126 


Figure 7- Harmonic response for fc = 0.05. 
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The results for k = 0.1 are shown in figure 8. The errors are quite small for this case with 
N pc = 320 and wave length A > 20. The total time T max = 20t r, the minimum allowed by the 
guidelines. 


MODEL EQUATION FOR CAP-TSD 

Quintic 32 pts, k = 0.100 

k ReF ImF 

0.000 - 0.0001 0.0000 
0.100 0.9983 -0.0045 

0.200 0.0003 0.0000 



I I I I I 


Implicit, 2nd order dt 

dt = 0.09817 k = 0.100 

Newton = 1 

0.0100 < dz < 1.0123 

0.0970 < CN < 9.8175 

Npc - 320 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 640 



4 0 z 20 



0 + 63 


Figure 8- Harmonic response for Jfe = 0.1. 
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The results for k = 0.25 are shown in figure 9. The time step Af = 0.09817 is still at its 
upper limit and N pc = 128, which is twice the minimum suggested. Five cycles of oscillation are 
required to meet the minimum T max requirement. More than one and one-half wave lengths are 
present across the mesh, which results in about 20 mesh points per wave length. 

The Fourier components are still quite accurate with the exception of the imaginary part of 
the first harmonic which is in error by about 1.5 percent. This term indicates a small phase error 
in the time history during the last cycle. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, k = 0.250 


k 

0.000 

0.250 

0.500 

0.750 

1.000 


ReF 

- 0.0002 


ImF 

0.0000 


1.0004 -0.0155 


-0.0008 

- 0.0012 

-0.0016 


0.0001 

0.0008 

0.0018 


Implicit, 2nd order dt 

dt = 0.09817 k = 0.250 

Newton = 1 

0.0100 < dz < 1.0123 

0.0970 < CN < 9.8175 

Npc = 128 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 640 


/A 


Figure 9- Harmonic response for k 



The results for k = 0.5 are shown in figure 10. For this case, N pc has attained its minimum 
value of 64. The dispersion error is revealed in the slower propagation of the numerical wave 
shown at the upper right. This slowdown leads to a significant difference between the numerical 
and exact spatial wave lengths. Here, there are approximately 6 mesh points per wave length in 
the center of the mesh, too few to maintain spatial accuracy. This inaccuracy leads to internal 
mesh reflections which result in an error of over 4 percent in the fundamental Fourier coefficent. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, k = 0.500 

k ReF ImF 

0.000 0.0013 0.0000 

0.500 0.9581 -0.0039 

1.000 -0.0073 0.0003 



J 


Implicit, 2nd order dt 

dt = 0.09817 k = 0.500 

Newton = 1 

0.0100 < dz < 1.0123 

0.0970 < CN < 9.8175 

Npc = 64 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 640 



Lun 1 I 1 1 1 l 1 JL-i 1 J- i II ,1 I Jill i il 

4 0 z 20 



I i i i f i i i i i i i i i i i i i i i i t i j i i i i i i i j i i i i i i i i I i i i i i i i i i i i i l i i i i i i i i i ] 

0 + 63 


Figure 10- Harmonic response for fc = 0.5. 
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The results for k = 0.75 are shown in figure 11. For this case At has been reduced to 
maintain N cy at its minimum value of 64. The mesh has approximately 8 points per wave length 

tk 2 7 *" d , 4 points near z = 10. These are the minimum values recommended above. 

Therefore, this frequency is the highest one for which the mesh fineness might be considered 
adequate. This relative mesh coarseness results in noticeable inaccuracy in the time history 

? f by dU ® 10 b0th internal mesh re,lections (e-9-. for t » 30) and boundary reflections (for 
• ^h 4 °l' ? er6 ' S a , s,gn,ficant error in the wave propagation speed across the mesh and the error 
in the fundamental Fourier coefficent now exceeds 8 percent. Although these errors are larqer 
than those which would be acceptable for this model problem, the corresponding errors for a TSD 

application would not be expected to be as large because of the reduction in amplitude which 
occurs for multidimensional wave propagation. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, k = 0.750 

k ReF imF 

0.000 - 0.0098 0.0000 

0.750 0.9543 - 0.0871 

1.500 0.0173 0.0170 



Implicit, 2nd order dt 

dt = 0.06545 k = 0.750 

Newton = 1 

0.0100 < dz < 1.0123 

0.0647 < CN < 6.5450 

Npc = 64 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 960 



mm i i i i i I 


0 


'I 

20 


o 



I- 1 1 1 1 1 1 1 i i i i i 

0 





J 


j i i i m 
63 


Figure 1 1- Harmonic response for k = 0.75. 
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The results for k — 1.0 are shown in figure 12. For this case, the mesh is obviously too 
coarse for the wave to propagate correctly and the spatial plot shows the typical saw-tooth wave 
with about 2 points per wave length which results. As was pointed out in the discussion of figure 4 
above, this mesh cannot propagate a wave of this frequency with any accuracy. The amplitude 
is much too small (at this instant of time when t = T max ). The time history shows large errors, 
including a phase error of about one-quarter cycle at r majt . This error is present even though the 
boundary condition is forcing the downwash to continue its sinusoidal oscillation. 


MODEL EQUATION FOR CAP-TSD 

Quintic 32 pts, k = 1.000 

k ReF ImF 

0.000 0.0517 0.0000 

1.000 -0.0306 -0.8457 

2.000 -0.0665 0.0432 



Implicit, 2nd order dt 

dt = 0.04909 k = 1.000 

Newton = 1 

0.0100 < dz < 1.0123 

0.0485 < CN < 4.9087 

Npc = 64 

Zmax = 20. Nz = 32 

Tmax = 63. Nmax = 1280 



O 

CL 



4 


turn t i i i t i t l i i t I 1 t I 1 1 1 1 1 lI 
0 z 20 



f 


! 



The results for k = 1.5 are shown in figure 13. For this high frequency case (i.e. p high for 
the mesh used), the wave cannot propagate across the mesh, the disturbances remain trapped 
near the origin, and the time history contains large apparently random errors. 


MODEL EQUATION FOR CAP-TSD 
Quintic 32 pts, k = 1.500 


0.000 

-0.1234 

0.0000 

Npc = 

64 

1.500 

0.3078 

-0.3962 

Zmax = 

20. 

3.000 

0.0494 

-0.0498 

Tmax = 

63. 

4.500 

0.0157 

-0.0481 

H A A 

A A 

6.000 

0.0084 

-0.0352 

■■■■ 

h A 


Implicit, 2nd order dt 

dt = 0.03272 k = 1.500 

Newton = 1 

0.0100 < dz < 1.0123 

0.0323 < CN < 3.2725 

Npc = 64 

Zmax = 20. Nz = 32 
Tmax = 63. Nmax = 1920 


Figure 13- Harmonic response for k = 1.5. 

In order to verify the usefulness of the guidelines for mesh parameter selection, a new mesh 
was designed to correct the complete breakdown in the solution for k = 1.5 seen in figure 13. 
The mesh design program of reference 2 was used with far-field and maximum spacings chosen 
to meet the criteria 

^far-field = w/8 k — 0.26 

/i max = 7r/4 k = 0.52 

for Jb = 1.5. This new mesh required N Zmxx = 61 points, nearly twice as many as the 32 points 
used previously. 






The results for Jb = 1.5 on the finer mesh are shown in figure 14. The maximum error in 
the fundamental Fourier coefficient is about 9 percent and the overall quality of the calculation is 
comparable to that obtained on the original mesh at k = 0.75, as seen in figure 1 1. As compared 
with the coarse mesh calculation in figure 13, the wave can now propagate across the mesh and 
the severe internal mesh reflections have been eliminated. Although not shown herein, the results 
on the 61 point mesh were improved over those on the 32 point mesh at all frequencies. 


MODEL EQUATION FOR CAP-TSD 

Quintic 61 pts, k = 1.500 

k ReF ImF 

0.000 0.0104 0.0000 

1.500 1.0210 0.0903 

3.000 - 0.0041 0.0036 



Implicit, 2nd order dt 

dt = 0,03272 k = 1,500 

Newton = 1 

0.0100 < dz < 0.5190 

0.0631 < CN < 3.2725 

Npc = 64 

Zmax = 20. Nz = 61 




I I I I t 1 1 I I i i 1 1 - 1 I I I i t 1 i I 1 I I I L 
0 



I I I I I 1 i 1 t 1 I I I t I I I 1 i 

63 


Figure 14.- Harmonic response with improved mesh for k = 1.5. 
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Although the maximum frequency k = 1.5 for which calculations were shown is perhaps 
much higher than would be of interest for flutter, it has been assumed that the Mach number is 
approximately 0.7 throughout the presentation of these results. An increase in the Mach number 
for the linear problem or significant transonic effects would reduce the frequency range for which 
the conclusions are valid. Said differently, a finer mesh would be needed to produce comparable 
results as the Mach number was increased. 

Conclusions 

Recommendations for choosing some of the input parameters used in the CAP-TSD computer 
code have been made. These parameters include those associated with the mesh design and time 
step. These guidelines are based primarily on experience with a one-dimensional model problem 
used to study wave propagation in the vertical direction. Numerical simulations presented herein 
support the recommendations given. 
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